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Abstract. - We study the decoherence of a superconducting qubit due to the dispersive coupling 
to a damped harmonic oscillator. We go beyond the weak qubit-oscillator coupling, which we 
associate with a phase Purcell effect, and enter into a strong coupling regime, with qualitatively 
different behavior of the dephasing rate. We identify and give a physicaly intuitive discussion 
of both decoherence mechanisms. Our results can be applied, with small adaptations, to a large 
variety of other physical systems, e. g. trapped ions and cavity QED, boosting theoretical and 
experimental decoherence studies. 



Introduction. — With a thrust from applications in 
quantum computing, the manipulation of quantum states 
in superconducting nanocircuits has made tremendous 
progress over the last decade [1-9] . A crucial step for 
these successes is the understanding of decoherence and 
the design of good measurement schemes. The latter is a 
particular challenge as the detector is made using the same 
technology as the system being detected i.e. the qubit. 
Also, the measurement timescale cannot be considered to 
be infinitesimally short as compared to the intrinsic scales 
of the qubit evolution. Thus, understanding the measure- 
ment process is crucial both fundamentally and for im- 
proving experiments. 

A specifically attractive development is the emergence of 
circuit quantum electrodynamics (cQED) [10-17], where 
effective Hamiltonians, similar to those of the coherent 
light-matter interaction of quantum optics and in partic- 
ular of cavity QED, can be realized in the microwave fre- 
quency domain. There are many approaches to realize the 
qubit, including flux and charge, and the cavity, including 
a superconducting quantum interference device (SQUID) 
or a coplanar waveguide. 

In this context, measurement protocols making use of 
dispersive qubit-oscillator interactions [1,2] are useful for 
reducing the backaction on the qubit [18]. For example. 



in the flux qubit-SQUID combination, as in the Delft 
setup of Refs. [1, 19], the SQUID behaves like a har- 
monic oscillator. Its inductive coupling to the flux qubit 
leads to a frequency shift depending on the qubit state 
ri|,l = vfP±A^. Here, Vl is the bare oscillator frequency 
and A is the quadratic frequency shift. A measurement 
of the SQUID resonance frequency provides information 
of the qubit state. While the manipulation of the qubit is 
usually performed at the optimum working point [3], the 
readout can and should be performed in quantum nonde- 
molition measurement i.e. in the pure dephasing limit. 

In this letter we study the decoherence of a qubit due to 
the dispersive coupling to a damped harmonic oscillator, 
taking the Delft setup as an example though our results 
may be adapted to several physical systems. In the Purcell 
effect a narrow oscillator linewidth enhances the absorb- 
tion of the resonant photon emitted by the two-level atom 
and thus the energy relaxation of the latter. In the weak 
qubit-oscillator couphng regime (WQOC), we explain the 
behavior of dephasing in terms of a similar process, the 
phase Purcell effect. This regime is characterized, as we 
wiU be show later, by A/f] < ^/njn / {I + n{n)Y / , where 
ri(ri) is the Bose function at the frequency Q, and envi- 
ronment temperature T. The main result of this work 
lays beyond the WQOC, in a regime where fast qubit- 
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oscillator entanglement plays the dominant role. We find 
a qualitatively diflcrcnt behavior of the dcphasing rate. 
The divergence of the qubit dephasing rate l/r^ oc 1/k 
when the oscillator decay rate k — > is lifted by the onset 
of the strong coupling regime. 

The Hamiltonian describing the Delft setup [19] can be 
written as 



Caldeira-Leggett way by a bath of harmonic oscillators 



H = 



E 



-Uz + ( a) a 



HA' 
4fl 



Hs 



(1) 



D 



Here, a and a) are the annihilation and creation operators 
of the harmonic oscillator, ct^ acts in the Hilbert space of 
the qubit and Hd describes the damping of the oscillator. 
A full-length derivation of Hamiltonian (1) and discus- 
sion of the approximations used is given in Ref. [20]. It 
basically derives the Hamiltonian from the equations of 
motion of the Josephson phases across the junctions and 
truncates the SQUID potential to the second order. 

We will show that key experiments [1, 2] are per- 
formed outside the WQOC. Moreover, a very recent ex- 
periment [21] explicitly relies on the use of a strong dis- 
persive coupling regime. We demonstrate that the dephas- 
ing rate l/r^ cx 1/k for WQOC, and I/t^ cx k at strong 
coupling. We discuss the crossover between these regimes 
and its dependence on k and temperature T. We provide 
physical interpretations of both regimes, the former as a 
phase Purcell effect and the latter as the onset of qubit- 
oscillator entanglement. The results of the present study 
may be extended straightforwardly to any system with 
similar dispersive qubit-oscillator coupling: the charge- 
qubit-coplanar wave guide system (see Yale setup [2]), 
trapped ions [22] and 3D microwave cavity QED [23], 
quantum dots [24] , among others. 



Method. — In studying the qubit dephasing we are 

facing the challenge of a complex non-markovian environ- 
ment consisting in the main oscillator (i.e. SQUID) and 
the ohmic bath. Moreover, the qubit couples to a non- 
Gaussian variable of its environment. Therefore the tools 
developed for Gaussian baths [25] cannot be applied in this 
system for arbitrary strong coupling between the qubit 
and the oscillator. 

We study the qubit dynamics under the Hamiltonian (1) 
for arbitrary A/$7, assuming essentially the dimensionless 
oscillator decay rate k/0 as the only small parameter. In 
this regime we avoid over-damping of the oscillator and the 
strong backaction on the system which this would cause. 
We give in the following a brief description of the crucial 
steps and approximations of the calculation. We model the 
damping, associated with the oscillator decay rate n, in the 
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with J{uj) = J^j ^'jfi/ ('2TnjU!j)S{uj — luj) — mhKUjQluj — 

uJc)/t^ and He the counter term [26-28] where 6 is the 
Heaviside step function and Wc an intrinsic high frequency 
cut-off. Our starting point is the Born-Markov master 
equation in the weak coupling to the bath limit for the 
reduced density matrix ps in the qubit-oscillator Hilbert 
space 
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This approach is valid at finite temperatures ksT ^ hn, 
for times t » l/wc [28,29], which is the limit we will dis- 
cuss henceforth. We start from a standard factorized ini- 
tial state for all subsystems. We express ps{t) in the qubit 
basis and represent its elements, which are still oscillator 
operators, in phase-space as 



PS 
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where Xo-o-' is the characteristic Wigner function and 
D = exp(aa^ — a* a) the displacement operator [30]. In- 
dependent of our work, Ref. [31] has used a different 
phase-space representation to calculate the qubit dephas- 
ing rate. We characterize the qubit coherence by C{t) = 
(o-a;(g)i) = 2Re Tr/5|^(t) which can be easily shown to 
be C{t) = 87rRex||(0, 0, t). After a rather long but es- 
sentially straightforward calculation, one obtains for xu 
a generalized Fokker-Planck equation 

X^l{a,a*,t) = ((^a{ki + in) + a* ki) da 
+ (a* {k2 - in) + ak2) da* 

-'-^{da - da')^ + p(a + a*f)xn(a,a*,^), (3) 

where 

fci,2 = -j(2T^(l + 2nT)±^(l + 2n,)), (4) 

^ = -^(%(l + 2nT)+J^t(l + 2ni))-^ (5) 

and Her = n{na) is the Bose function. To solve Eq. (3) we 
make a Gaussian ansatz for xti- 

XTi = A{t) exp{-M{t)a^ - N{t)a*^ - Q{t)aa*).{6) 
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This ansatz includes coherent and thermal states. In the 

following wc assiimc the oscillator to be initially in a ther- 
mal state, in equilibrium with its environment. This 
implies Q(0) = 1/2 + n{n) and M(0) = N{0) = 0. Due 
to the quadratic (pure dcphasing) form of the Hamilto- 
nian (1), obtain a closed system of ordinary differential 
equations for the parameters of the Gaussian ansatz, see 
also Ref. [20] . This system can be easily solved perturba- 
tively in A in the weak coupling regime, or numerically, 
(for arbitrarily strong coupling), and we can extract the 
dephasing time from the strictly exponential long-time 
tail of C(t) = 8nRcA{t). 

Weak qubit-oscillator coupling. — Before solving 
Eq. (3) in a general manner, we revisit the case of small A. 
Up to the lowest non- vanishing order A'^, the analytically 
calculated WQOC dephasing rate is 



= A' 



4n(0) (n(0) + l) 



02 



(7) 



where k„ = ^2kBTn/{h{l + 2n{n))). The term 1/k 
exactly reproduces the Golden Rule dephasing rate of 
Ref. [19], and is similar to the result of Ref. [32]. These 
previous results have been obtained considering only the 
two-point correlator of the fluctuating observable {a+a'^)'^, 
i.e. assuming an Gaussian environment. The crossover 
point Km from 1/k to k in Eq. (7) is, at the Delft param- 
eters [1], comparable to O, i.e., k would dominate over 
1/k only in a regime where the Born approximation fails. 
Nevertheless, since the golden rule limit limK-»oo ^/t<I> = 
is unphysical, such a term was to be expected. 

In the WQOC regime, the enhancement of dephasing 
by weak coupling to the environment is analogous to the 
enhancement of spontaneous emission by the narrow cav- 
ity lines in the resonant Purcell effect, see Refs. [33,34]. 
In the pure dephasing case we have no energy exchange 
between the qubit and the oscillator. Qubit decoherence 
is caused by fluctuations of (a + a^)^. Since we are in the 
WQOC regime, the stronger coupling between the oscilla- 
tor and the environment causes equilibrium between the 
oscillator and the bath on a shorter time scale than the 
qubit dephasing. In equilibrium, the main contribution to 
the fluctuations of {a + a^)"^ is the exchange of photons 
between oscillator and bath. The process is analogous to 
equilibrium fluctuations in canonical thermodynamics. A 
virtual photon returning from the environment is at res- 
onance with the oscillator. The absorbtion of this pho- 
ton, like in the resonant Purcell effect will be enhanced 
by narrow oscillator lines. Therefore, the entire dephasing 
process will be enhanced when the coupling to the envi- 
ronment is weak and this mechanism can be viewed as a 
phase Purcell effect. We give a more detailed discussion 
of this effect in the Appendix. 

Strong qubit-oscillator coupling. — The dephas- 
ing rate (7) obtained in the small k and WQOC limit 
diverges for k 0, i.e., in the absence of an environment. 



The solution to this apparent contradiction lies beyond 
the WQOC, therefore wc solve Eq. (3) numerically using 
again the Gaussian ansatz for xu- 
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Fig. 1: Dephasing rate l/r^ as function of A for different val- 
ues of K. Power-law growth at low A crosses over to A- 
independence at strong coupling. Inset: Dephasing rate as a 
function of k in the weak coupling regime (A/n = IQ-*) show- 
ing l/r^ oc A*/k and the strong coupling regime {A/fl = 10~^) 
with 1/t^ oc k dependence. Here hQ/kBT = 2 similar to ex- 
periments. 

Figure 1 shows the dependence of the dephasing rate on, 
A for various values of k. The dimensionless parameter 
h^l/ksT is 2, similar to the Delft and Yale setups. As pre- 
dicted by eq. (7) for k <^ Km and small A, the dephasing 
rate is proportional to A^/k. Further increasing A, we 
observe a saturation of the dephasing rate which marks 
the onset of the strong coupling regime. This regime is 
analogous to the strong coupling in linear cavity QED. 
Here 1 /r^ is proportional to k. 

At strong qubit-oscillator coupling the oscillator couples 
to the qubit stronger than it couples to the heat bath, such 
that one cannot use the effective bath concept of WQOC. 
As the qubit-oscillator system becomes entangled, a fun- 
damentally different dephasing mechanism emerges. The 
eigenstates of the Hamiltonian (1) at k = are the dressed 
states {)(T, m^)} where \mcr) are the number states of the 
oscillator with frequency fio-- Opposed to WQOC, these 
dressed states are built in the strong coupling regime on a 
shorter time scale than the re-thermalization of the oscilla- 
tor. In the evolution from thermal state of oscillator with 
frequency f2 to an equilibrium between the new oscillator 
with and the bath, the state in the narrower potential 
tends to absorb and the one in the wider potential to emit 
photons to the bath in an incoherent manner, causing fluc- 
tuations of (a + a^)^ and thus qubit decoherence. Thus we 
expect 1/t^ oc Kn{Q). This simple picture is confirmed 
by numerical results in Fig. 2, for a wide range of values 
of K. The inset of Fig. 2 shows the crossover from strong 
coupling rate k to WQOC rate 1/k. This indicates that, 
for fixed A, as k decreases, A stops being "small" and the 
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Fig. 2: Scaling plot of the dephasing rate I/t^ as function of 
temperature. For A/f2 = 0.3, i.e. in the strong coupling regime 
(see Fig. 1), for a wide range of k's wc show that 1/{t^k) is 
proportional to the Bose function n(Q,). Inset; dephasing rate 
l/r^ as function of k for different values of A and HQ/kBT = 2. 
Continuous lines correspond to Kn(f2), dashed lines correspond 
to n(n)(n(n) + l)A4/(n2/t). 

WQOC limit breaks down. Thus, approaching k = for 
any given A we eventually leave the domain of validity for 
eq. (7) avoiding the divergence at k — > 0. As expected, 
dephasing will vanish as wc go to a finite quantum system 
(qubit (g) single oscillator) at k = 0. We observe that the 
criterion of " small" A in WQOC is valid only relative to k. 
Using 1/t0 = Kn{Q) in the strong coupling regime and the 
1/k term of I/t^ in eq. (7) in the weak coupling regime, 
we determine the position of the crossover Ac between the 
two regimes 



A, 



1 



f2(l + n(f2))V4 



(8) 



The position of the cross-over is controlled by the ratio 
of the coupling strengths between the three subsystems 
i.e. A^/f2 and k. Note that, with the in-situ tuning of 
the qubit-SQUID coupling, available in the Delft experi- 
ment, the position of the cross-over could be tested exper- 
imentally. Using the parameters from Ref. [1,2] one finds 
(A/Ac)Yaic ~ 1-4 and (A/Ac)Dcirt ~ 1-3 i.e., the strong 
coupling regime finds application in both setups . 

If the oscillator is weakly driven off-resonance, as is the 
case in the dispersive measurement, the qualitative behav- 
ior remains the same as in Fig. 1, as shown in Ref. [20]. 
In general a tunneling term may occur in Eq. (1) and 
lead to energy relaxation as well as further reducing the 
matrix elements containing the dephasing rate. We ex- 
pect that, as long as the energy splitting E of the qubit is 
off-resonance with the oscillator, which in our case means 
\E — 2n\ ^ K, the effect of the relaxation is rather weak 
and dephasing still dominates. On resonance, we expect 
a similar Purccll to strong coupling crossover as for the 
dephasing channel. 



Our results have applications in other systems with 
similar dispersive qubit-oscillator coupling, e.g., the Yale 
setup [2], in the off-resonant dispersive regime. There, the 
system is described by a similar (Eq. (12) in Ref. [32]) 
quach'atic; coupling a^a between qubit and cavity and a 
pure dephasing Hamiltonian. In particular, a strong dis- 
persive regime of this system has been utilized to resolve 
number states of the electromagnetic field in Ref. [21]. 
The terms and a'^'^ in Eq. (1) do not play a cen- 
tral role for our physical predictions, as confirmed by the 
numerical simulations. Wc expect our results, with mi- 
nor adaptations, to be applicable to various cavity sys- 
tems, e.g. quantum dot or atom-based quantum optical 
schemes [23,24]. The dispersive coupling of Hamiltonian 
(1) could have implications for the generation of squeezed 
states, quantum memory in the frame of quantum infor- 
mation processing, measurement and post-selection of the 
number states of the cavity. 

Conclusion. — We have presented a concise theory of 
the dephasing of a qubit coupled to a dispersive detector 
spanning both strong and weak coupling. The phase-space 
method applied is based on treating the oscillator on the 
same level of accuracy as the qubit. We have discussed 
the dominating decohercncc mechanism at weak qubit- 
oscillator coupling, where the linewidth of the damped 
oscillator plays the main role, analogous to the Purcell ef- 
fect. At strong qubit-oscillator coupling we have identified 
a qualitatively different behavior of the qubit dephasing 
and discussed it in terms of the onset of the qubit-oscillator 
entanglement. We have provided a criterion delimitating 
the parameter range at which these processes dominate 
the qubit dephasing. 
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Appendix. — Assuming the WQOC limit we use 
Fermi's golden rule in an otherwise exact manner to prove 
the analogy between the weak qubit-oscillartor coupling 
regime and the Purcell effect. One can map the damped 
oscillator by an exact normal mode transformation [35] 
onto an effective heat bath of decoupled oscillators denoted 



by Cj , & and with a spectral density 



Jeff(a;) = 



2ku) 



(^2 _ Q2)2 + ^2^^ 



2,, ,2 ■ 



(9) 



Jeff corresponds to the effective density of electromagnetic 
modes in the cavity introduced in regular linear cavity 
QED for describing the Purcell effect. The WQOC deco- 
herence rate is proportional to the two-point correlation 
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function of the environmental operator coupling to the 
qubit [28,36], in our case 

S2{^) = {x\t)X'iO))^ - ^{X'))\ (10) 

where X is the sum of the effective bath coordinates 
X = J2j -\/?i/(2mjWj)(cj + Cj). For the pure dephasing 
situation described by the Hamiltonian (1) we only need 
to study l/r^ oc S2{co — > 0+) because the qubit energy 
conservation implies energy conservation within its effec- 
tive environment. The last term of Eq. (10) removes the 
noise bias. This is important since dephasing is caused 
only by processes that leave a trace in the bath [37], i.e. the 
exchanged boson spends a finite time in the environment. 
Terms of the structure ^ct(t)ct(t)cfeC;^ ,(^Ci{t)cj{t)cl.cl'^ 
contribute to 5*2(0) only when LUi = ujj = 0, which arc 
modes with density J^g — 2ku)/^'^ each, leading to terms 
are of order k^. Up to linear order in k, the only terms in 
S2{u) — > 0+) that fulfill the energy conservation and leave 

a trace in the bath are of the structure (|c|(t)cj(t)c|c;^, 

including the permutations among the operators taken at 
time t and those taken at time 0. The terms contribut- 
ing to 52 (w 0+) satisfy the condition \uji — ujj\ 0+. 
Physically this corresponds to infinitesimal energy fiuctu- 
ations which leave a trace in the bath. Or, in other words, 
the photon absorbed at t=0, cj, should spend finite time 
in the bath and be emitted back only at the later time 
t, but at the same time the energy change in the envi- 
ronment e.g. caused by cjc; should remain undetectable 
within the energy-time uncertainty at every time, there- 
fore in the Golden Rule (long time) limit wi « u)j . Taking 
the continuum limit we thus have 

/•OO 

1/t0 cx / Jeff (u;)(l + «(tj)) Jeff (u;)n(w). (11) 
Jo 

The integral in eq. (11) can be rewritten as the convolution 
K{oj') = J (hjJes{uj)n{u!)Jes{'^' — oj)n{u)' — oj) for co' 0. 
Using eq. (9), K{ui') becomes a function with resonances 
at oj' = and uj' = 2^1. The the height of these resonances 
and consequently l/r^ oc 52(0) increases with decreasing 
K, thus matching the behavior of the dephasing rate (7). 
At the same time, the tail of the peak at 2f7 enhances 
52(0) when K increases. This corresponds to the k term in 
eq. (7). Analogous to l/r^ in eq. 7, 52 (w — > 0+) vanishes 
for T ^ 0. 
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